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Abstract: We study the non-perturbative properties of = 2 super conformal field 
theories in four dimensions using localization techniques. In particular we consider SU(2) 
gauge theories, deformed by a generic e-background, with four fundamental flavors or with 
one adjoint hypermultiplet. In both cases we explicitly compute the first few instanton 
corrections to the partition function and the prepotential using Nekrasov's approach. These 
results allow to reconstruct exact expressions involving quasi-modular functions of the bare 
gauge coupling constant and to show that the prepotential terms satisfy a modular anomaly 
equation that takes the form of a recursion relation with an explicitly e-dependent term. 
We then investigate the implications of this recursion relation on the modular properties of 
the effective theory and find that with a suitable redefinition of the prepotential and of the 
effective coupling it is possible, at least up to the third order in the deformation parameters, 
to cast the S-duality relations in the same form as they appear in the Seiberg-Witten 
solution of the undeformed theory. 
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1. Introduction 

Four-dimcnsional field theories with rigid M = 2 supersymmetry provide a remarkable 
arena where many exact results can be obtained; indeed M = 2 supersymmetry, not being 
maximal, allows for a great deal of flexibility but, at the same time, is large enough to guar- 
antee full control. This fact was exploited in the seminal papers [1, 2] where it was shown 
that the effective dynamics of = 2 super Yang-Mills (SYM) theories in the limit of low 
energy and momenta can be exactly encoded in the so-called Seiberg-Wittcn (SW) curve 
describing the geometry of the moduli space of the SYM vacua. When the gauge group 
is SU(2), the SW curve defines a torus whose complex structure parameter is identified 
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with the (complexified) gauge coupUng constant r of the SYM theory at low energy. This 
coupling receives perturbative corrections at 1-loop and non-perturbative corrections due 
to instantons, and the corresponding effective action follows from a prepotential T that is 
a holomorphic function of the vacuum expectation value a of the adjoint vector multiplet, 
of the flavor masses, if any, and of the dynamically generated scale in asymptotically free 
theories or of the bare gauge coupling constant tq in conformal models (see for instance [3] 
for a review and extensions of this approach). 

Recently, M = 2 superconformal field theories (SCFT) have attracted a lot of atten- 
tion. Two canonical examples of SCFT's are the A/" = 2 SU(2) SYM theory with iV/ = 4 
fundamental hypermultiplets and the J\f = 2* theory, namely a J\f = 2 SYM theory with 
an adjoint hypermultiplet. In both cases, the /3-function vanishes but, when the hypermul- 
tiplets are massive, the bare coupling tq gets renormalized at 1-loop by terms proportional 
to the mass parameters. Besides these, there are also non-perturbative corrections due to 
instantons. As shown in [4] for the J\f = 2* theory and more recently in [5] for the Nj = 4 
theory, by organizing the effective prepotential ^ as a series in inverse powers of a and by 
exploiting a recursion relation hidden in the SW curve, it is possible to write the various 
terms of J- as exact functions of the bare coupling. These functions are polynomials in 
Eisenstein series and Jacobi ^-functions of tq and their modular properties allow to show 
that the effective theory at low energy inherits the Sl(2, Z) symmetry of the microscopic 
theory at high energy. In particular one can show [4, 5] that the S-duality map on the 
bare coupling, i.e. tq — > — I/tq, implies the corresponding map on the effective coupling, 
i.e. r ^ — 1/r, and that the prepotential T and its S-dual are related to each other by a 
Legendre transformation. 

The non-perturbative corrections predicted by the SW solution can also be obtained 
directly via multi-instanton calculus and the use of localization techniques [6, 7]^. This ap- 
proach is based on the calculation of the instanton partition function after introducing two 
deformation parameters, ei and e2, of mass dimension 1 which break the four-dimensional 
Lorentz invariance, regularize the space-time volume and fully localize the integrals over 
the instanton moduli space on sets of isolated points, thus allowing their explicit evaluation. 
This method, which has been extensively applied to many models (see for instance [10] - 
[17]) can be interpreted as the effect of putting the gauge theory in a curved background, 
known as Jl-background [6, 7, 18], or in a supergravity background with a non-trivial 
graviphoton field strength, which are equivalent on the instanton moduli space [19, 20]. 
The resulting instanton partition function ^inst(ei, eg), also known as Nekrasov partition 
function, allows to obtain the non-perturbative part of the SYM prepotential according to 

^inst = - lim eie2 logZinst(ei,e2) . (1.1) 

Actually, the Nekrasov partition function is useful not only when the e parameters are 
sent to zero, but also when they are kept at finite values. In this case, in fact, the non- 
perturbative e-deformed prepotential 

-^mst(ei,e2) = -eie2 log Zinst(ei, £2) (1-2) 
^See also [8, 9] for earlier applications of these techniques. 
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represents a very interesting generalization of the SYM one. By adding to it the corre- 
sponding (e-deformed) perturbative part -Fpert ) one gets a generahzed prepotential that can 
be conveniently expanded as follows 

oo 

i^pert + Finst = ^ ^^"''^ + ^2)'" (^1^2)^ • (1-3) 

n,g=0 

The amplitude F^^''^^ which is the only one that remains when the e-dcformations are 
switched off, coincides with the SYM prepotential of the SW theory, up to the classical 
tree-level term. The amplitudes with g > 1 account instead for gravitational cou- 

plings and correspond to F-terms in the effective action of the form F^^'^^W'^^ , where W 
is the chiral Weyl superfield containing the graviphoton field strength as its lowest com- 
ponent. These terms were obtained long ago from the genus g partition function of the 
M = 2 topological string on an appropriate Calabi-Yau background [21] and were shown 
to satisfy a holomorphic anomaly equation [22, 23] which allows to recursively reconstruct 
the higher genus contributions from the lower genus ones (see for instance [24, 25]). More 
recently, also the amplitudes with n 7^ have been related to the J\f = 2 topolog- 

ical string and have been shown to correspond to higher dimensional F-tcrms of the type 

is a chiral projection of real functions of = 2 vector superfields 
[26], which also satisfy an extended holomorphic anomaly equation [27, 28]. By taking 
the limit €2 — )• with ei finite, one selects in (1.3) the amplitudes F^^'^\ This limit, also 
known as Nekrasov-Shatashvili limit [29], is particularly interesting since it is believed that 
the M = 2 effective theory can be described in this case by certain quantum integrable 
systems. Furthermore, in this Nekrasov-Shatashvili limit, using saddle point methods it is 
possible to derive a generalized SW curve [30, 31] and extend the above-mentioned results 
for the SYM theories to the e-deformed ones. 

By considering the Nekrasov partition function and the corresponding generalized pre- 
potential for rank one SCFT's, in [32] a very remarkable relation has been uncovered with 
the correlation functions of a two-dimensional Liouville theory with an e-dependent central 
charge. In particular, for the SU(2) theory with Nj = 4 the generalized prepotential turns 
out to be related to the logarithm of the conformal blocks of four Liouville operators on 
a sphere and the bare gauge coupling constant to the cross-ratio of the punctures where 
the four operators are located; for the M = 2* SU(2) theory, instead, the correspondence 
works with the one-point conformal blocks on a torus whose complex structure parameter 
plays the role of the complexified bare gauge coupling. Since the conformal blocks of the 
Liouville theory have well-defined properties under modular transformations, it is natural 
to explore modularity also on the four-dimensional gauge theory and try to connect it to 
the strong/weak-coupling S-duality, generalizing in this way the SW results when ei and 
£2 are non-zero. On the other hand, one expects that the deformed gauge theory should 
somehow inherit the duality properties of the Type IIB string theory in which it can be 
embedded. Some important steps towards this goal have been made in [33] and also in 
[34] where it has been shown that the SW contour integral techniques remain valid also 
when both ei and €2 are non-vanishing. To make further progress and gain a more quanti- 
tative understanding, it would be useful to know the various amplitudes F^"'^) in (1.3) as 
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exact functions of the gauge coupling constant and analyze their behavior under modular 
transformations, similarly to what has been done for the SW prepotential F^^'^^ in [4, 5]. 
Recently, by exploiting the generalized holomorphic anomaly equation, an exact expression 
in terms of Eisenstein series has been given for the first few amplitudes F^^'S) of the SU(2) 
theory with Nf = 4 and the J\f = 2* SU(2) theory in the limit of vanishing hypermul- 
tiplet masses [35]. This analysis has then been extended in [36] to the massive M = 2* 
model in the Nekrasov-Shatashvili limit, using again the extended holomorphic anomaly 
equation, and in [37] using the properties of the Liouville toroidal conformal blocks in the 
semi-classical limit of infinite central charge. However, finding the modular properties of 
the deformed prepotential in full generality still remains an open issue. 

In this paper wc address this problem and extend the previous results by adopting a 
different strategy. In Section 2, using localization techniques we explicitly compute the 
first few instanton corrections to the prepotential for the N = 2* massive theory with 
gauge group SU(2) in a generic e-background. From these explicit results we then infer 
the exact expressions of the various prepotential coefficients and write them in terms of 
Eisenstein series of the bare coupling. Our results reduce to those of [4, 5] when the 
deformation parameters are switched off, and to those of [35] - [37] in the masslcss or in 
the Nekrasov-Shatashvili limits. The properties of the Eisenstein series allow to analyze 
the behavior of the various prepotential terms under modular transformations and also 
to write a recursion relation that is equivalent to the holomorphic anomaly equation if 
one trades modularity for holomorphicity. The recursion relation we find contains a term 
proportional to 6162, which is invisible in the SYM limit or in the Nekrasov-Shatashvili 
limit. In Section 3 we repeat the same steps for the M = 2 SU(2) theory with Nf = 4 
and arbitrary mass parameters and also in this case derive the modular anomaly equation 
in the form of a recursion relation. In Section 4 we study in detail the properties of the 
generalized prepotential under S-duality, and show that when both ei and 62 are different 
from zero, due to the new term in the recursion relation, the prepotential and its S-dual are 
not any more related by a Legendre transformation, an observation which has been recently 
put forward in [38] from a different perspective. We also propose how the relation between 
the prepotential and its S-dual has to be modified, by computing the first corrections in 
eie2. Finally, in Section 5 we conclude by showing that there exist suitable redefinitions 
of the prepotential and of the effective coupling that allow to recover the simple Legendre 
relation and write the S-duality relations in the same form as in the undeformed SW theory. 
The appendices contain some technical details and present several explicit formulas which 
are useful for the computations described in the main text. 

2. The N = 2* SU(2) theory 

The M = 2* SYM theory describes the interactions of a A/" = 2 gauge vector multiplet 
with a massive M = 2 hypermultiplet in the adjoint representation of the gauge group. 
It can be regarded as a massive deformation of the = 4 SYM theory in which the 
/?-function remains vanishing but the gauge coupling constant receives both pcrturbative 
and non-perturbative corrections proportional to the hypermultiplet mass. Using the lo- 
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calization techniques [6, 7] one can obtain a generalization of this theory by considering the 
e-dependent terms in the Nekrasov partition function. In the following we only discuss the 
case in which the gauge group is SU(2) (broken down to U(l) by the vacuum expectation 
value a of the adjoint scalar of the gauge vector multiplet). We begin by considering the 
non-perturbative corrections. 

2.1 Instanton partition functions 

The partition function at instanton number k is defined by the following integral over 
the instanton moduli space Mk- 



= J dA1fee"'^'°=* (2.1) 



where 5inst is the instanton moduli action of the M = 2* theory. After introducing two 
deformation parameters ei and €2, the partition function Z]^ can be explicitly computed 
using the localization techniques. In the case at hand, each Z^ can be expressed as a sum 
of terms in one-to-one correspondence with an ordered pair of Young tableaux of V{k) 
such that the total number of boxes in the two tableaux is k. For example, at /c = 1 
we have the two possibilities: (□,•) and (•,□); at A: = 2 we have instead the five cases: 
(□□,•), (•,□□), (0)*)) (•jB)) (□)□); and so on and so forth. Referring for example to 
the Appendix A of [39] for details, at A; = 1 one finds 



(2.2) 



^ __ (-ei + m)(-e2 + m){au + m){a2i +m — ei - €2) 

^■^'•^ (-ei)(-e2)ai2(a2i - ei - 62) 

^ _ (— ei -I- m)(-e2 + m){a2i + 'm){ai2 + m — ei — €2) 

^'''-'^ (-ei)(-e2)a2i(ai2 - ei - €2) 

where auv = o,u — O'v with ai = —02 = a, and 

~ ei + £2 o\ 
m = m-\ ^ — (^--jj 

is the equivariant mass parameter for the adjoint hypermultiplet in the e-background [40, 
35]. The 1-instanton partition function is therefore 

7 _ 7 _ (W - (ei - e2)^)(16a^ - 4m^ - 3(6i + €2?) ,^ 
^1 - ^(n,.) + 8eie2(4a2-(ei + e2)2) ' ^^'^^ 

At A; = 2 the relevant partition functions are 



Z, 



_ (—61 -I- m)(— €2 + m)(-ei + €2 + m)(-2e2 + m)(ai2 + m) 
" (-ei)(-e2)(-ei + e2)(-262)ai2 

^ (ai2 + m + £2) (021 + m - ei - £2) (021 + m - ei - 2e2) 
(ai2 + e2)(a2i - ei - e2)(a2i - ei - 2e2) 

_ (-ei + m)^(-e2 + w,)^(ai2 + fh - ei)(ai2 + fh- €2) 
" (-ei)2(-e2)2(ai2-ei)(ai2-e2) 
^ (021 -f m - ei)(a2i + m - 62) 
(a2i - ei)(a2i - €2) 



(2.5) 
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The contributions corresponding to the other Young tableaux at k = 2 can be obtained 
from the previous expressions with suitable redefinitions. In particular, is obtained 

from ■^(□□,,) in (2.5) by exchanging ai2 -H- a2i; ) obtained by exchanging ei -H- €2, 

and finally Z r— 1 is obtained by simultaneously exchanging ai2 -H- 021 and €1 €2- The 

complete 2-instanton partition function is 

^2 = -^(m,*) + -^(•,m) + ^(Q,.) + + -^(□,n) (2-^) 

but we refrain from writing its expression since it is not particularly inspiring. This pro- 
cedure can be systematically extended to higher instanton numbers leading to explicit 
formulas for the instanton partition functions. 

Following [7] we can cast these results in a nice and compact form. Indeed, defining 

q = e'^'^ (2.7) 

where tq = ^ + i ^ is the complexified gauge coupling constant of the J\f = 2* SYM 
theory, the grand-canonical instanton partition function 

oo 

Zi^,t = ^q^''Zk (2.8) 



fe=o 



where Zq = 1, can be rewritten as 

oo 2 

^mst = E n n 

{Yi,Y2) i,j=lu,v=l 



a-uv + ei(i - 1) - e2j 



auv + ei(i - 1 - kyj) - e2(j - kui) 
auv + m + - I - Kj) - e2(j 



X 



auv + m + ei(i - 1) - €2] 



(2.9) 



where the first line represents the contribution of the gauge vector multiplet and the second 
line that of the adjoint hypermultiplet. Here /c^j and k^i denote, respectively, the number 
of boxes in the i-th row and in the i-th column of a Young tableaux 1^ and are related to 
the Dynkin indices of the corresponding representation. These quantities can be extended 
for any integer i with the convention that kui = or = if the i-ik tow or the i-th 
column of is empty. For example, for Yy, = UJ the kui's are (2, 0, 0, 0, . . .) and the fc„i's 
are (1, 1, 0, 0, . . .). Moreover we have 

u,i u,i 

It is quite straightforward to check that the expressions (2.2) and (2.5) are reproduced by 
the compact formula (2.9) by selecting the appropriate Young tableaux. 

Following Nekrasov's prescription, we can obtain the generalized non-perturbative pre- 
potential according to 

00 

i^inst = -£162 log ^inst = ^ ^'^^ ■ (2-11) 

fe=l 
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In the limit — t- 0, the above expression computes the instanton contributions to the 
prepotential of the = 2* SYM theory, while the finite e-dependent terms represent further 
non-perturbative corrections. Notice that in the limit m — >■ 0, (2.9) correctly reduces to 
the partition function of the J\f = 4 SYM theory, namely to the Euler characteristic of the 
instanton moduli space (see for instance [41]). 

2.2 Perturbative part 

The compact expression (2.9) allows to "guess" the perturbative part of the partition 
function by applying the same formal reasoning of Section 3.10 of [6]. Indeed, in (2.9) we 
recognize the following universal {i.e. ^-independent but a-dependent) factor 

TT n '^nv + eiji - 1) - €2j 

which, if suitably interpreted and regularized [6, 7], can be related to the perturbative part 
of the partition function of the M = 2* theory in the e-background. According to this idea, 
we are then led to write 

-^pert = eie2 > > log -r. rr : . (2.13 

auv + m + ei{i - 1) - €23 

U,V=1 l,J = l \ / ^ 

Using the following representation for the logarithm 

X d / dt 



where A is an arbitrary scale and summing over i and j, we can rewrite (2.13) as 

2 



(2.14) 



u,v=\ 



-^pert = eie2 ^ \lei,e2{auv) " lei,e2{auv + m) (2.15) 

where (see also [18, 35]) 
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(2.16) 



'■^''^ ~ TsyT{s)Jo T(e-«*-l)(e-2*-l)J 

This function, which is related to the logarithm of the Barnes double F-function, can be 
easily computed by expanding for small values of ei and €2- As a result, i^pert becomes a 
series in inverse powers of whose first few terms are 

Fpe. = ^(4m^ - S^) log ^ - ^ L 

(4m^ - s^) (4m^ + 4p - s^) (4m^ + 6p - 3s^) 

61440 

(4m2 - s^) (4m2 + 4p - s^) (48m^ + I76m'^p + 160^^ - SSm^s'^ - 204ps'^ + 51s^) 

8257536 

+ 0{a-^). (2.17) 
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Here we have used (2.3) and introduced the convenient notation 



s = ei + 62 , p = 6162 . (2-18) 

One can easily check that in the hmit — t- this expression correctly reproduces the 1-loop 
prepotential of the M = 2* SU(2) gauge theory (see for instance [4]). 

2.3 Generalized prepotential 

The complete generalized prepotential is the sum of the classical, perturbative and non- 
perturbative parts. We now focus on the latter two terms which are directly related to 
the Nekrasov partition function. Just like the 1-loop piece (2.17), also the instanton terms 
(2.11) can be organized as a series expansion for small values of ei and £2, or equivalently 
for large values of a. Discarding a-independent terms, which are not relevant for the gauge 
theory dynamics, we write 

2 °° h \ 

Fpert + Fmst = ^^0 log ^ - ^ -^^ (2.19) 

where the coefficients he, are polynomials in m^, and p which can be explicitly derived 
from (2.17) as far as the perturbative part is concerned, and from the instanton partition 
functions (2.9), after using (2.11), for the non-perturbative part. Here we list the first few 
of these coefficients up to three instantons: 

/lo = ^(W -s2) ^ (2.20) 

h, = {4m' - s') {W + 4p- s') + ^ (2-21) 



h2 = (4m^ - s^) [Am^ + 4p - s^) 



3840 

_^ 12m' + 18p - 21s' ^4 ^ 8m' + 12p - 9s' 



+ + (2.22) 



16 

, o ON/. 2 . o W48m'^ + 176m' p+160p' -88m' s' -204ps' + 51s^ 
hs = {4m' - s') {4m' + 4p - s") ^ 

_ 3s^ 2 _ 240m4 + 1392m'p + 1440p2 - SOOOm^s^ - 7548^5^ + 3903s'' 4 
~"32"^ 1024 ^ 



240m^ + 1008m'p + 960^^ - 1080m's' - 2652ps' + 987s^ . 

32 ^ 



■) • (2-23) 



It is interesting to notice that all h^s are proportional to {4m' — s') and, except for hg, 
also to {4m' + 4p — s'). Exphcit expressions for h^ with £ > 3 can be systematically derived 
from the generalized prepotential but they are not needed for our considerations. 

Building on previous results obtained in — >■ limit from the SW curve [4, 5] and 
on the analysis of [35] for the massless case, we expect that the expressions (2.20)-(2.23) 
are just the first terms in the instanton expansion of (quasi) modular functions of q. More 
precisely, we expect that h^ are (quasi) modular functions of weight 2£ that can be written 
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solely in terms of the Eisenstein series E2, E^i and Eq (see Appendix A for our conventions 
and definitions). This is indeed what happens. In fact we have 

/io = ^(W-s'), (2.24) 

hi = (4m2 - s^) {Aw? + Ap - s^) E2 , (2.25) 

X [5(4m2 + 6p- s^)E^ + [Arri^ + 6p- 13s^)E^ , (2.26) 

= 214 . 33 . 5 . 7 (4"^' - (4™' + 4p - s') 

X [35(80m^ + 272m'^p + 2AQp^ - A{)rr?s^ - 68ps^ + hs'^)El 

+ 84(l6m^ + Urn^p + 6O/ - 56m^s^ - 136ps^ + 13s'')£;4 E2 

+ (I76m^ + MAm^p + 960^^ - ISlGm^s^ - 4556^5^ + 3323s'^)£;6j . (2.27) 

By using the small q expansion of the Eisenstein series one can check that the explicit 
instanton contributions we have computed using localization techniques are correctly re- 
covered from the previous formulas. We stress that the fact that the various instanton 
terms nicely combine into expressions involving only the Eisenstein series is not obvious 
a priori and is a very strong a posteriori test on the numerical coefficients appearing in 
(2.20)-(2.23). It is quite remarkable that the explicit instanton results at low k can be 
nicely extrapolated and allow to reconstruct modular forms from which, by expanding in 
powers of q, one can obtain the contributions at any instanton number. 

We can also organize the generalized prepotential according to (1.3) and obtain the 
amplitudes F^"'^) as a series in inverse powers of with coefficients that are polynomials 
in E2, E4 and Eq. The first few of such amplitudes are: 

^(0,0) _ 2 1 ^ _ ^^^2 _ m'^ibE^+E^) _ m^jUbEl + 84^2^4 + 11£;6) 

-m og^ 48 a2 5760 2903040 '"' 



(2.28) 



1 2a m^E2 m^iE^ + E^) {175E^ + 336E2Ei + 8Ee) 

"4 °^ a" ^ 96 a2 ^ 1536 a4 ^ 2903040 ^ 



^(0,1) ^ m'^E2 _ m*{J>El + E^) _ m^(llE| + 6E2F4 + ^e) 
' ~ 48 a2 2304 41472 



(2.29) 

+ • • • (2.30) 



p(2,o) ^ E2 _ m\bEl + 9E4) m\nbEl + 588ii;2i^4 + 559£;6) 

768 a2 30720 a4 7741440 ^ ' > 

f(^A) = m2(25ii;2 + 17£;4) m^(385Ff + 798i;2 + 213£;6) 

192 a2 23040 a4 1935360 ^ ' ^ 

(0,2) ^ _rrf{5El + E^ _ m^lGOEj + 93E2E4 + UE^) 

3840 414720 a6 ^ ' ' 
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Terms with higher values of n and g can be systematically generated without any difficulty 
from the /i^ 's given in (2.24)-(2.27). The first term F^^''^^ represents the prepotential of the 
M = 2* SU(2) gauge theory and its expression (2.28) agrees with that found in [4] from 
the SW curve (see also [5]). The other terms are generalizations of those considered in [35] 
and more recently in [37] in particular limits (m — > or 62 ^ 0) where they drastically 
simplify. 

2.4 Recursion relations 

The generalized prepotential (2.19) is clearly holomorphic by construction but does not have 
nice transformation properties under the modular group since the coefficients Tn explicitly 
depend on the second Eisenstein series E2 which is not a good modular function. Indeed, 
under 

TO ^ , with ad - 6c = 1 , (2.34) 

cTo + a 

E2 transforms inhomogeneously as follows 

6c 

E2{tq) ^ (cTo + dfE2{TQ) + (cro + d) — . (2.35) 



TTl 



Therefore, in order to have good modular properties we should replace everywhere E2 
with the shifted Eisenstein series E2 = E2 + at the price, however, of loosing 

holomorphicity. This fact leads to the so-called holomorphic anomaly equation [22, 23] 
(see also [28]). On the other hand, in the limit tq — s- 00, holding r fixed so that E2 E2, 
we recover holomorphicity but loose good modular properties and obtain the so-called 
modular anomaly equation [42]. This equation can be rephrased in terms of a recursion 
relation satisfied by the he coefficients which allows to completely fix their dependence on 

To this aim let us consider the explicit expressions (2.24)-(2.27) and compute the 
derivatives of /i^ with respect to E2 . With simple algebra we find 

|| = ^E/>A-.-. + ^.../.« (2.36) 

with the initial condition 

We have explicitly checked this relation for several values of i > 3; we can thus regard 
it as a distinctive property of the e-deformed N = 2* low-energy theory. Notice that in 
the limit — )• 0, (2.36) reduces to the recursion relation satisfied by the coefficients of the 
prepotential of the M =2* SU(2) theory found in [4, 5] from the SW curve, and that the 
linear term in the right hand side disappears in the so-called Nekrasov-Shatashvili limit 
[29] where one of the two deformation parameters vanishes and a generalized SW curve 
can be introduced [30, 31]. 

Eq. (2.36) can be formulated also as a recursion relation for the amplitudes F*^"'^) 
defined in (1.3). To see this, let us first extract the a-dependence and, in analogy with 
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(2.19), write 

^^"''^ = -E^^' (2-38) 



j:{n,g) ^ 



e=i 
SO that 



he = ^fr\ei + e2f^{eie2r . (2.39) 

n,9 

Notice that the coefficients which are polynomials in the hypermultiplet mass m, 

have mass dimensions 2(1 + £ — n — g); therefore f^"''^^ = Oiin + g>£ + l (this condition 
can be easily checked on the explicit expressions (2.28)-(2.31). These definitions must be 
supplemented by the "initial conditions" 

/r=m^ /r = -^, /r)=0 (2.40) 
which are obtained from (2.24). Inserting (2.39) in the recursion relation (2.36), we obtain 

^/i"'^^ _ ^ f(^^'9^) An2,92)\ , - 1) f{n,g-l) ... 

nl,n2;pi,p2 i=0 

where the ' means that the sum is performed over all ni, n2, gi and g2 such that ni+n2 = n 
and gi + g2 = g- Eq- (2.41) shows that the coefficients /^^"'^•' and hence the amplitudes 
jp{i^,9) recursively related to those with lower values of n and g, similarly to what can 
be deduced from the holomorphic anomaly equation [35] . 

3. The Af^2 SU(2) theory with A^/ = 4 

We now consider the Af = 2 SU(2) theory with Nf flavor hypermultiplets in the funda- 
mental representation of the gauge group. When Nf = 4 the 1-loop /3-function vanishes 
and the conformal invariance is broken only by the flavor masses ruf (/ = 1,...,4). Fur- 
thermore there are non-perturbative effects due to instantons which are nicely encoded in 
the exact SW solution [2] . We now discuss the generalizations of these effects induced by 
the €£ deformation parameters in the Nekrasov partition function, following the same path 
described in the previous section for the M = 2* theory. 

3.1 Instanton partition functions 

Using localization techniques, we can express the instanton partition functions of the 
Nf = 4 theory as sums of terms related to pairs of Young tableaux. Referring again the 
Appendix A of [39] for details, at = 1 we have the following two contributions 

1 ^ 
^(□,.) = (-e,)(-e.)a,.(a.r-er-e.) J^}^' + ' 

^ ; (3-1) 



(-ei)(-e2)a2i(ai2 - ei - £2) jj-^ 
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where 



mf = irif + 



ei + €2 



(3.2) 



is the equivariant hypermuhiplct mass in the e-background^. Summing the two terms (3.1) 

and setting ai = — 02 = o, we get 



■^1 - ^{n,») + ^{;n) 
1 



2eie2 



+ ^ m f m ji + (ei + €2) ^ + - (ei + €2)^ + 



4Pfm 



/</' 



/ 



4a2 - (ei + £2)2 



(3.3) 

where Pfm = 771177127713777,4. The partition function Zi contains a part proportional to a?, 
a part which does not depend on a and a part containing in the denominator. In all 
a-dependent terms the flavor masses 777/ always occur in S0(8)-invariant combinations and 
thus in all such terms we can always express the mass dependence using the following 
quadratic, quartic and sextic SO (8) invariants: 



R 



777) , 



/</' / 



(3.4) 



N = 



_3_ 
16 



777^777^,777^,, - — ^ 777^^777^, + ^ 2^ "7^ . 

/</'</" f^f f 

The a-independent terms in (3.3), instead, are not invariant under the 80(8) flavor group. 
Indeed, neither "^fiTif nor '^f^f/'n^fm'f' respect this invariancc. Notice, however, that 
o-independent terms in the partition functions, and hence in the prepotcntial, are irrele- 
vant for the gauge theory dynamics, and thus can be neglected. We will always do so in 
presenting our results. 

The explicit expressions of the partition functions Zf^ for k > 1 arc rather cumbersome; 
nevertheless it is possible to write the (grand-canonical) instanton partition function in a 
quite compact way using the connection with the Young tableaux. Indeed, denoting as x 
the instanton counting parameter and using the same notations introduced in Section 2.1, 
we have 



-'inst 



k=0 



E n 



{Yi,Y2) 



a-uv + ei(7 - 1) - e2j 



J^i o-uv + ei(« - 1 - Kj) - e2(j - Ki) 



(3.5) 



nn 

«=i/=i 



au + mf + 61 (7 - 1) - e2(j - Ki) 
au + rhf + ei(7 - 1) - e2j 



^Note that in [39] the hypermultiplet masses were denoted — m/. 
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Here the second line represents the contribution of the gauge vector multiplet and the 
last line that of the fundamental hypermultiplets. From this expression, by selecting the 
appropriate Young tableaux, one can obtain the various terms of the instanton partition 
function and their dependence on the parameters. 



3.2 Perturbative part 

Also in the Nf = 4 theory one can deduce the perturbative contribution to the partition 
function from the "universal" factor of the grand-canonical instanton partition function 
(3.5), namely 



n 



1 



2 2 4 

n,'"- ^ "'^ - - n n ^ _ ^ _ 



(3.6) 



which, of course, requires a suitable interpretation and regularization. Using this obser- 
vation and following again [6, 7], we can write the perturbative part of the prepotential 
as 



-^pert = eie2 



au + nif + ei{i - 1) - e2j 



u,v=l 



u=l f=l 



which, after using (2.14) and summing over i and j, becomes 

-^^pert = £162 



2 2 4 

u,v=l u=l f=l 



(3.7) 



(3. 



where the function j^-^^^ is defined in (2.16). Expanding the latter for small values of ei 
and 62 and discarding as usual a-independent terms, we obtain 



Fpert = - log 16 a' + ^ {4R -s^+p) log J 



1 



96 a2 
1 



3840 



16R^ - 96Ti - 8R{s^ - 2p) + - 4s^p + 3p^ 

128R^ - 1920i?ri + 768iV - 160{R^ - 6Ti){s'^ - 2p) 



(3.9) 



+ i?(56s^ - 224s^p + 192p2) - 6s^ + 36s^p - 63sV + 30p^ + O [a'^) 

where s and p are as in (2.18). One can easily check that in the limit — ^ this expression 
correctly reproduces the 1-loop prepotential of the M = 2 SU(2) Nf = 4 theory (see for 
instance [3]). Notice the presence of a term proportional to a? which corresponds to a finite 
1-loop renormalization of the gauge coupling constant. 



13 



3.3 Generalized prepotential 



We now consider the instanton corrections to the prepotential of the Nf = 4 theory which 
is given by 

oo 

i^inst = -6162 log ^inst = ^ Ffc . (3.10) 

k=l 

At = 1, up to a-independent terms, from (3.1) we simply get 

a2 2(ri + 2T2) 
F, - -eie2 Z,--- _ ,2 

The presence of an a^-tcrm in Fi signals that there is a rcnormalization of the gauge 
coupling constant at the 1-instanton level. Actually, the same thing occurs also at higher 
instanton numbers; indeed we find 

ISa^ ^ 23a2 2701 

F2 = h • • • , -F3 = H , -F4 = + • • • (3.12) 

^ 64 ' ^ 192 ' 32678 ^ ^ 

Combining these contributions with the perturbative one given by the first term in (3.8) 
and adding also the classical tree- level prepotential = log xa^, we obtain 

, 1 13 2 23 , 2701 4 \ 2 , 2 

i^\ogx-\ogl6 + -x + — x^ + —x^ + ^27gg + • • • ja' = logga^ (3.13) 

that is a redefinition of the instanton expansion parameter and hence of the gauge coupling 
constant. By inverting this relation we get [43] 

x = 16q(l-8q + Uq^ - 192q^ + •••) = tj^t^ (3-14) 

(see Appendix A for some properties of the Jacobi ^-functions). As pointed out for instance 
in [44], by setting 

q = e'^^^o with ro = - + i ^ , (3.15) 

one can show that tq is precisely the modular parameter of the SW curve for the SU(2) Nf = 
4 theory appearing in the original paper [2]. This gauge coupling constant receives only 
corrections proportional to the hypermultiplet masses and as such is the strict analogue^ of 
the gauge coupling constant tq of the M = 2* theory considered in Section 2. On the other 
hand, using (3.14) one can show that x is a cross-ratio of the four roots of the original SW 
curve [2]. Prom now on we always present the results for the Nf = 4 theory in terms of q. 



Except for a customary overall factor of 2. 
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Expanding the complete prepotential -Fpert + -^inst as in (2.19), the first coefficients 
up to three instantons turn out to be 



ho 
hi 



1 

24 



(3.16) 



16(i?2 - GTi) - 8i2(s2 - 2p) + 8^^- As^p + 3^^ + 32(Ti + 2T2) q 



h2 = ^ 



16(i?2 + GTi) - 8i?(s2 - 2p) + s'' - 4s2p + 3^^ 
+ 128(ri + 2r2) q^ + Oiq"^) , 

\28R^ - 1920iiri + 768iV - 160{R^ - 6Ti){s^ - 2p) 
+ i?(56s'^ - 224s2p + I92p2) _ ggS _^ 35^4^ _ 63^2^2 _^ 3q^3 
+ 32s2(ri + 2r2)g 



(3.17) 



+ 



768iV + 384ii:ri - 64(i?2 + e,Ti){s'^ -p) + m{4:S^ - Us'^p + 5p2) 



4s^ + 20s^j9 - 25s2p^ + 



128(Ti + 2T2) {IQR - lls^ + IQp) q^ + 0{q^ 



(3.18) 



We worked out also the expressions for a few other /i^'s with £ > 3 which however become 
rapidly very cumbersome and thus we do not write them here (see for example Appendix 
B where we give the detailed form of /;,,3). 

As in the N = 2* model, also in the Nf = 4 theory we can view the previous results as 
the first instanton terms of the expansion of (quasi) modular forms in powers of q. In this 
case, however, besides the Eisenstein series £'2, £'4 and Eq, also the Jacobi ^-functions 62 
and 04 are needed. Matching the g-expansion of these modular functions with the explicit 
results (3.16)-(3.18), guided by what we already obtained in the — )• limit [5] which 
suggests the existence of a recursion relation that fixes the E2 dependence, we find^ 



ho=\{m-s^+p) , 

hi = i^(4i2 - s'^+p) {4R -s^ + 3p) E2 - 4(ri - T2 9^) 

h2 



24 

1 



(3.19) 
(3.20) 



— (4i? -s^+p) {AR -s^ + 3p) [AR - + Ap) El 



- {AR - s2 + 4p) (Ti ej - T2 9l) E2 



+ 



720 



6AR^ - 80R^{3s^ - Ap) + 4ii(27s'' 



- 13s^ + 68s^p - Ms^p^ + 30/ + 2304Ar 



8s> + 49^"= 
£4 



^[R-s^+p) [ti ei{2ei + ei) + T2 et{ei + 2ei) 



(3.21) 



*As compared to the results presented in [5, 44], here all masses have been rescaled by a factor of \/2, 



that is mf"'" = ^/2i 



for all /. 
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We have checked that a similar pattern occurs also in other /i^'s with £ > 3 (see Appendix 
B). 

Organizing the complete prepotential as an expansion in powers of and p as in (1.3), 
we can obtain the amplitudes F^"'^) , the first few of which are 



Fm = 9RW^ - ^ ^ T,9i-T2ei _ R\hEl+E,) _ NE, 
6a2 ^ a2 180 Sa^ 

RT,ei{2E2 + 2el + ei) RT20t{2E2 + 2ei + ei) 

+ 6a4 6a4 ^"^"^^^ 

i.(i,o) - + R\El + E,) 

2^A^12a2^ 48a4 

_ T^ei{E2 + Aei + 2et) T2el{E2 - Aoj - 2el) 

12 + 12 a4 ^^"^"^^ 

f(o,i) = i ,o " _ _ ^'(2i^l + £^4) 

2 ^ A 6a2 36a4 

^ rig|(2g2 + 291 + g|) r2g|(2E2 - 2ei - el) ^ 

6 6 ' 



^(2,0) ^ _ _ ^^(5^1+9^4) 



96a2 960 a4 ' ^^'^^^ 



= _^ + i?(10i^| + 11^4) 



24^' 360^^^ + ^'-''^ 

pio,2)_ E, R(95Ei + 1<)E,) 

~ 32a2 2880 a4 + ' 

pm) _ 5EI + 13^4 j.(2,i)_ ioi^i + i7£;4 . . 

11520 ' 2880 ' ^ ' 

^(,,2)^95^+94^... j,(o.3)_ '^El + E, 

11520 ' 384 ■ ^ ' ^ 

One can easily check that F^'*^) in (3.22) completely agrees with the prepotential of the 
M = 2 SU(2) Nf = 4 gauge theory as derived for example in [5]; notice also that in the 
massless limit, i.e. R,N,Ti,T2 — >■ 0, these amplitude drastically simplify and precisely 
match the results presented in [35]. It is also interesting to observe that in the Nckrasov- 
Shatshvili limit where one of the e^'s vanishes, the amplitudes F*^"'*^-* of the Nj = 4 theory 
reduce to those of the J\f = 2* theory upon setting Ti = T2 = = and R = m? . Indeed, 
with these positions and rescaling a — )■ 2a, the amplitudes F^'*''^) in (3.22)-(3.29) become 
twice the corresponding amplitudes of the M = 2* theory given in (2.28)- (2.33). This 
simple relation does not hold away from the Nekrasov-Shatashvili limit: the amplitudes 
F^"''?) with g 7^ are in fact intrinsically different in the two theories, as a consequence of 
their different gravitational structure. 



3.4 Recursion relations 

Looking at the explicit expressions (3.19) - (3.21), it is not difficult to realize that the E2 
dependence of the h^s is quite simple and exhibits a recursive pattern. This points again to 
the existence of a recursion relation among the h^s involving their derivatives with respect 
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to E2. Indeed, following the same steps described in Section 2.4, one can check that 



= Q hiht-i-1 H €162 (3.30) 

with the initial condition 

This recursion relation has exactly the same structure of that of the M = 2* theory given 
in (2.36), the only difference being in the numerical coefficients. Notice that the coefficients 
of the quadratic terms can be matched by a rescaling of the /i^'s, but those of the linear 
terms proportional to 6162 remain different for the two theories; this is another signal of 
their intrinsically different behavior when a generic e-background is considered. 

Finally, if we expand the h^s as in (2.39) we can reformulate the recursion relation 
(3.30) in terms of the partial amplitudes /f"'^^ and get 



Qf{n,9) ^ , ^-1 

ni,n2;5i,g2 j=0 



n ^ 

dE2 6 



E' ( E /i'^"'^ ViT!^ + ftV"-^ (3.32) 



where the ' means that the sum is performed over all ni, n2, g\ and 52 such that ni+n2 = n 
and gi + g2 = g- This relation has to be supplemented by the initial conditions 

/r=2i?, fr'=-i, fr'=\ (3.33) 

obtained from (3.19). 

In the next section we will analyze the implications of the recursion relations and in 
particular their consequences on the modular transformation properties of the generalized 
prepotential. 



4. Modular anomaly equations and S-duality 

The fact that the generalized prepotential can be written in terms of (quasi) modular 

functions of the bare coupling constant allows to explore the modularity properties of 
the deformed theory and study how the Sl(2, Z) symmetry of the microscopic high-energy 
theory is realized in the effective low-energy theory. In the following we will concentrate on 
the SU(2) Nf = 4 theory, even if our analysis and our methods can be equally well applied 
to the J\f = 2* SU(2) model. Furthermore, we will focus on the generator S of the modular 
group, corresponding to the following transformation of the bare coupling constant 

S : TO ^ - - . (4.1) 

To 

Notice that this action implies that the instanton counting parameter x in the Nekrasov's 
partition of the Nf = 4 theory, defined in (3.14), transforms as follows 

5 : X ^ 1-x , (4.2) 
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as one can readily check from the properties of the Jacobi ^-functions. This type of trans- 
formation is consistent with the interpretation of x as a cross-ratio, a fact that is also 
exploited in the AGT correspondence with the Liouville conformal blocks [32]. As dis- 
cussed in [2], in the Nf = 4 theory the modular group acts with triality transformations 
on the mass invariants (3.4); in particular one has 

S : R ^ R , Ti ^ T2 , T2 ^ Ti , N ^ N . (4.3) 

In the deformed theory, these rules have to be supplemented by those that specify how 
Sl(2, Z) acts on the equivariant deformation parameters. Adapting to the present case the 
considerations made in [45] for the e-dcformcd conformal Chern-Simons theory in three 
dimensions, we assume that S simply exchanges ei and 62 with each other^. In particular 
this means that s = ei + €2 and p = €162 are left unchanged, i.e. 

S : s ^ s , p ^ p . (4.4) 

Using the rules (4.1) - (4.4) and the modular properties of the Eisenstein series and Jacobi 
^-functions, it is easy to show that the coefficients hi of the generalized prepotential (see 
(3.21) and Appendix B) transform as quasi-modular forms of weight 2£ with anomalous 
terms due to the presence of the second Eisenstein series E2, namely 

S : k,XE,)^ rf ft, + J-) = g i (^) ' D>'HAE.) (4.5) 

where we have introduced the convenient notation D = 129^2 • These transformation rules 
are formally identical to those of the underformed Nf = 4 theory derived in [5] from the 
SW curve; however, since the hg's satisfy a modified recursion relation with a new term 
proportional to eie2, the practical effects of (4.5) in the deformed theory are different from 
those of the undeformed case, as we shall see momentarily. 

To proceed let us consider the pair made by a and its S-dual image S{a) = an, on 
which S acts as follows [1, 2]: 






S: L 1 ^ r. ; 1 L I = ( - I . (4.6) 
We therefore have 

S^{a) = -a . (4.7) 
In the SW theory this relation is enforced by taking 

5(a) = (4.8) 
^ ' 2tt\ da ^ ' 



®Note that this rule is consistent with the interpretation of ei and €2 as fluxes of (complex) combinations 
of NS-NS and R-R 3-form fleld strengths in a Type IIB string theory realization which rotate among 
themselves under S-duality. 
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where T is the undeformed effective prepotential which is related to its S-dual by a Legendre 
transform: 

T -S{T) = 1'K\aS{a) . (4.9) 

It seems natural to try the same thing also in the deformed theory. Recalling that the 
effective generalized prepotential is 

oo , ^ 

F = TTiro a' + holog J- J2^,-ri (4-10) 

e=i 

where the first term is the classical part, we therefore posit 

Applying the S-duality rules (4.1) and (4.5), we obtain 

Sia) 1 f^f^e^j^^h^^O^ 



1 

-a + 



' Dhi I ^ {2i + l)hehr^ 



(27riTo) 



2^ a'^W Z-^ 2^+" Q2£+2n+3 

e=o e,n=o 



The expression in square brackets can be simplified using the recursion relation (3.30); like 
in the undeformed theory, the quadratic terms in the /I's exactly cancel but, due to the new 
e-dependcnt term in the deformed recursion relation, an uncanceled part proportional to 
€162 remains. This simple calculation shows that in order to enforce the relation (4.7) when 
€162 7^ 0, the standard definition (4.11) has to be modified by adding terms proportional 
to €i€2 in the right hand side. In the sequel we will work oTit explicitly the first corrections 
and show how the relation (4.7) constrains the form of S{a). 

4.1 S-duality at first order in £162 

To organize the calculation, we introduce a set of generating functions ipe for the coefficients 
h£ according to 

^0 = -K\og- + ^^^:^—^ , ipi+i = -daipi for e>0. (4.13) 
e=i " 

In particular we have the following relations with the generalized prepotential F: 

(po = Triro - F , (4.14a) 
(pi = -27riro a + daF , (4.14b) 
cp2 = 27riro - d^F = 27riro - 27rir , (4.14c) 

where in the last line we have introduced the effective coupling r. 
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As shown in Appendix C, the ip^s form a ring under the action of D due to the gen- 
erahzed modular anomaly equations (3.30). This ring structure will enable us to formally 
express everything as functions of the c^^'s with coefficients that may depend on the product 
€162. Such a dependence is a consequence of the eie2-term in the recursion relation (3.30) 
and is the only explicit dependence on the deformation parameters that will be relevant 
for our purposes, all the rest being implicit inside the (^^'s and the hiis therein. 

As argued above, the definition (4.11) for S{a), namely S{a) = tqu + </?i/(27ri), has to 
be replaced by a new one containing terms proportional to €162, i.e. 



with 



oia) = Tn a H 

^ ^ 27ri 27ri 



X = eie2Xi + (£162)^X2 + 



(4.15) 



(4.16) 



where the X^'s have to be determined by imposing the constraint (4.7). Applying S-duality 
to (4.15), it is straightforward to obtain 



-a + 



1 



27riro - 



{roSiipi) - ipi) + {toS{X) - X) 



(4.17) 



therefore, in order to satisfy the relation (4.7), the expression in the square brackets above 
must vanish. Expanding this condition in €162, we obtain the following constraints 



(ro5(<^i) - (pi) 



, 



(ro S{^i) - (^1) + V {to S{Xk) - Xk) 



k=l 



n—k 



(4.18a) 
(4.18b) 



where the symbol |^ means taking the coefficient of (6162)" (notice that the (/p^'s defined in 
(4.13) do not have any explicit dependence on 6162 and thus ip^\^ = 5nQ for the same 
reason wc also have |^ = bno X^. 

Let us now compute 5((^i). Using the S-duality rules (4.1) and (4.5), and exhibiting 
temporarily the dependencies on E2 and a which are the only relevant ones for our purposes, 
we obtain 



ro5(</?i(£^2;a)) = ^ 



^=0 



= "Pi ( 



^2 + 



2^ (5(a)) 
6 



2^+1 



\ 2^+1 



-; a + 



+ — 
27riro 27ri/ 



(4.19) 



TTlTo 

= e^^-°+^^"Vi(-E2;a) 
Actually, this is a particular case of the more general result 



r^<S((^^) = e^-i-o 



(4.20) 
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Expanding the exponential and using the relation 



n 



(4.21) 



after some straightforward algebra to rearrange the various terms, we obtain 

-1)" / X \n. 



n=0 



n! 



V 2771X0/ 



(4.22) 



where the functions S^^^ are defined by 



C=¥'i 



As shown in Appendix C, these functions satisfy the remarkably simple relation 



(4.23) 



(4.24) 



that is a consequence of the ring properties obeyed by the functions ip^, which in turn are 
due to the modular anomaly equations (3.30). In view of this we can therefore rewrite 
(4.22) as 

(4.25) 



ro^5(<^,)=e2fc^»EW 



Note that the bare coupling tq initially appearing in the right hand side of (4.20) has been 
dressed into the effective coupling r. 

We now exploit this result and proceed perturbatively in 6162 to obtain explicit ex- 
pressions. At the zeroth order in 6162, from (4.25) we have 



T0<S((^1 



(4.26) 



where the last equality follows from (C.13). The constraint (4.18a) is therefore identically 
satisfied. This is no surprise since we have already shown that the relation (4.7) must be 
true up to terms proportional to 6162. At the first order in the deformation parameters we 
get instead 



tqS{(Pi) 



1 27rir 



Xi ip2 

~ 2(27ri)r2 27rir 



(4.27) 



where we have used (C.13). Inserting this result into (4.18b) for n = 1, we obtain the 
following equation for Xi 





S{Xi 



(4.28) 



r 2(27ri)r2 ' 

To solve it we take advantage of the S-duality relations at the zeroth order in 6162, which 
are formally the same of the SW theory. Using (4.24) - (4.26) it is easy to show that 



TqT 







(4.29) 
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and also that 

Sir) 1 = --. (4.30) 
Equipped with these results, one can check that a solution to (4.28) is given by 

^' = &- <^-^^ 

In conclusion^ we find that the S-dual image of a, which obeys the constraint (4.7) up to 
terms of order (£162)^, is 

5(a) = ro a + |i + e,e, + OHe.e^f) • (4.32) 

Observing that (^3 = liri daT, we can rewrite the eie2-term above also as 

'^\9,logf-) , (4.33) 



4(27ri) Vro/ 
so that (4.32) becomes 

2i7r5(a) = daF+'-^dalog (^^^ +0{{eie2f) . (4.34) 

In the following we will investigate the implications of this result, while we refer to Ap- 
pendix D for its extension to the second order in eie2- 

4.2 S-duality on the prepotential 

An immediate consequence of the e-correction in (4.34) is that S{a) is not simply pro- 
portional to the derivative of the prepotential F; thus it is natural to expect that in the 
deformed theory the S-dual of the prepotential is not simply given by a Lcgcndre transfor- 
mation as it happens instead in the undeformed SW case (see (4.9)). In [38] the relation 
between F and S{F) has been conjectured to be a deformed Fourier transformation. Here 
we reach the same conclusion, even though from a different perspective since for us all 
deviations from the undeformed theory are parametrized by the product eie2. 

We can compute the first e-corrcctions to the relation between F and S{F) using the 
same methods of the previous subsection and the results given in Appendix C. The starting 
point is the relation between ipQ and the generalized prepotential (see (4.14a)); from this 
we get the useful identity 

TTl 9 

F-S{F) = 2-KiaS{a) + —iS{a) - aro) + S{ipo) - ipQ . (4.35) 

To 

Using (4.15) and (4.22) for £ = 0, we can easily rewrite the right hand side and obtain 

F - 5(F) = 2„ „ + ^ + - ^0 - (S"> - ^^^^ 



®Note that actually the solution (4.31) is not unique, since in principle one could add to Xi a term Y\_ 
such that S{Y\) = Y\/t. A detailed analysis shows that a term of this type is of the form afzlij^ -\- 1) 
with a constant, which has different pole structure in r with respect to X\. We therefore do not consider 
this possibility and take a = 0. 
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where the dots stand for terms of order X'^ which are at least C'((eie2)'^). Notice that 
the difference T,^^^ — (pi is of order €i€2, as one can see from (4.26); thus the hnear term 
in X gives contributions C'((eie2)^), Hke the term. This means that the knowledge 
of the first order correction Xi obtained in the previous subsection is enough to compute 
the correction to F — S{F) at order (eie^)'^. This mechanism actually works at all orders, 
namely the k-th order coefficient of F — S{F) does not depend on Xk but only on Xj with 
j < k. 

Using the explicit value of Xi given in (4.31) and the expressions for the S'^^-^'s given 
in Appendix C (see in particular (C.12)- (C.14)), after straightforward algebra we obtain 

F - SiF-, = 2.i„5(.) - '-f log {L) + + + . 

(4.37) 

This shows that the simple Legendre transform relation (4.9), which holds in the under- 
formed theory, does not work any more when £162 ^ 0, as argued, from a different point of 
view, in [38]. In the next section, however, we will show that with a suitable redefinition 
of the prepotential F and of the coupling constant r it is possible to recover the standard 
Legendre transform relation also in the eie2-deformed theory. 



5. Conclusions 

The explicit first-order calculation of Section 4.1 shows that S{a) is not simply the deriva- 
tive of the prepotential F. However, it is still a total derivative, as is clear from (4.34). 
This feature is maintained also at the second order. Indeed, as shown in Appendix D, the 
(6162)^ correction X2 can be chosen as 



1 (y95 23 <f3(p4 11 (fil 



16 (27rir)2 96 (27rir)3 64 (27rir)4 

= 8 (- - — \ 

"V 16(27rir)2 192 {2mTf) ' 

so that we can rewrite (4.15) in the following way 



(5.1) 



S{a) = 2^ daF (5.2) 

with 

In the deformed theory the S-duality transformation of the effective coupling r does not 
have a simple form; in fact, using (4.14c) and the transformation properties of (p2, it is 
easy to show that (see (D.7)) 
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However, there is a modified effective coupling on which S-duality acts in a simple way. 
This is 



r = daS{a) 



27ri " 



eie2 



ip4 



+ 



4 V(27rir)2 (27riT 



(5.5) 



One can easily check that 



1 

r 



(5.6) 



Thus, in the effective e-deformed theory it is r, and not r, that exhibits the same behavior 
of the bare coupling tq under S-duality. Wc also observe that if one uses r, the expression 
of the extended prepotential F given in (5.3) simplifies and becomes 



i^ + ^log 



'eie2\2 



To) ( 4 ) 12(2 



TTir 



(5.7) 



This result seems to suggest that it is possible to write the higher order corrections only 
in terms of ips which, being the triple derivative of the prepotential, is proportional to the 
Yukawa coupling Caaai the rank-three symmetric tensor playing a crucial role in special 
geometry. Moreover, one can verify the simple Legendre transform relation 



F-S{F) = 2TTiaS{a) 



(5.8) 



up to terms of order (€162)^. It is quite natural to expect that this pattern extends also to 
higher orders. 

Our detailed analysis shows that when both deformation parameters ei and €2 are 
non-vanishing, besides the e-dependent structures generated by the Nekrasov partition 
functions, the effective theory seems to require a new series of explicit eie2-corrections in 
order to have S-duality acting in the proper way. These new corrections, being proportional 
to inverse powers of the coupling constant, appear to correspond to perturbative terms at 
higher loops and are absent in the Nekrasov-Shatashvili limit. We find remarkable that 
by using the modified prepotential F and the modified coupling f , all S-duality relations 
acquire the standard simple form as in the undeformed theory (see (5.2), (5.6) and (5.8)). 
It would be very nice to see whether these results admit an interpretation in the context 
of special geometry or in a more general geometric set-up that allows to go beyond the 
perturbative approach in the deformation parameters we have used in this paper. It would 
be interesting also to study the recursion relations obeyed by the functions ip/s which can 
be iteratively obtained from the viscous Burgers equation satisfied by ipi (see (C.2)). 

We conclude by observing that all S-duality formulas we have derived for the SU(2) 
Nf = 4 theory can be obtained for the J\f = 2* SU(2) theory as well. The only difference in 
this case is that every explicit occurrence of £162 has to be replaced by (eie2)/2 to take into 
account the different normalization of the viscous term in the modular anomaly recursion 
relation (see (2.36) as compared to (3.30)). 
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A. Moduleir functions 

We collect here some useful formulas involving the modular functions we used, 
^-functions: The Jacobi ^-functions are defined as 

e[l]{v\T) = ^g("-f)%2.Un-i)(.-|) ^ (A.l) 

for a,b = 0,1 and q = e^'"^. We simplify the notation by writing, as usual, 6i = 0[^], 
02 = 9^^] , 63 = 6 [q] , ^4 = ^ [^] . The functions 9a, a = 2,3, 4, satisfy the ^^aequatio identica 
satis abstrusa" 

0l _0l_0j = o^ (A.2) 
and admit the following series expansions 

92{0\t) = 2qi{l+q'^ + q'^ + q^^ + ■■■) , 

93{0\t) = l + 2q + 2q^ + 2q^ + --- , (A.3) 
94{0\t) = 1 - 2g + - + • • • . 

ry-function: The Dedekind ry-function is defined by 

00 

ri{q)=q^^ll(l-q'n ■ (A.4) 
n=l 

Eisenstein series: The first Eisenstein series can be expressed as follows: 

00 

E2 = l- 24^ai(n)g2" = 1 _ 2Aq'^ - 72g^ - 96^^ + ... , 

n=l 

00 

E4 = 1 + 240Y^ 0-3 (n) = 1 + 240g^ + 2160g^ + 6720g^ + . . . , (A.5) 

n=l 

00 

E6 = l-b04j2 ^5(") g^" = 1 - 504g2 - 16632^'^ - 122976g^ + . . . , 

n=l 

where ak{n) is the sum of the k-th power of the divisors of n, i.e., ak{n) = X^dln'^'^- '^^^ 
series E4 and Eq are related to the ^-functions in the following way 

E4 = l {9l + 9l + 9l) , E, = ^ {9l + 9i) {Ol + 61) {Oj - Ol) . (A.6) 

The series E2, E4 and Eq are connected among themselves by logarithmic g-derivatives 
according to 

qdgE2 = ^ {El - E4) , qdgE4 = ^ [E4E2 - Eq) , qdgEe = EqE2 - Ej . (A. 7) 



25 



Also the derivatives of the functions 9^ have simple expressions: 

qd.et = I (£^2 + 91 + et) , qd.el = | (£;2 + el - ej) , qd.et = | (^2 - ei - et) . 

(A.8) 

Modular transformations: Under the Sl(2, Z) modular transformation 
, ar + 6 



T ^ T 



with a,b,c,d E Z and ad — be = 1 , 



(A.9) 



CT + d 

the Eisenstein series -E4 and Eq behave as modular forms of weight 4 and 6, respectively: 
E4{t') = {cT + df E4{t) , Eeir') = {cr + df Eq{t) . (A. 10) 



The series E2, instead, is a quasi modular form of degree 2: 

E2{t') = (cr + df E2{t) + ^cicT + d) . 

in 



(A.ll) 



Under the generators T and S of the modular group the ^-functions and the Dedekind 77 
function transform as follows 



r : 

S : 



'4 ! 



ei2 ?7 



(A.12) 



el 



-r (^4 , 



^ 1^2 5 



B. The coefficient hs of the SU(2) A^/ = 4 prepotential 

Here we give the explicit expression of the coefficient up to three instantons obtained 
using localization. It is given by 



h3 



1 



1344 



768i?^ - 21504i?2ri + 12288i?Ar + 26112r^ - GUATiT^ - GUAT^ 



- (1792i2^ - 26880i?Ti + 10752N){s^ - 2p) + R^{1568s^ - 6272s'^p + 5376^^) 

- Ti(9408s^ - 37632s2p + 32256^^) - i?(496s^ - 2976sV + 5248sV - 2560^^) 

+ 51/ - 408sV + 1102sV - 1144sV + 357/1 + 24s'^(ri + 2X2)5 



3840(Ti + 2T2f - 1152(i?Ti + 2N){5s^ - 2p) 



+ 48(i?^ + 6ri)(16s* - 22s^p + 7p^) - 192R{2s^ - Ss'^p + 7s V - p ) 
+ 3(16s^ - 96s^p + 168sV - 88sV + 9/)j q'^ 
+ 32(ri + 2T2) [l60(i?2 + 6T1) - 2R{280s'^ - 272p) 
+ SlSs'' - 856s2p + 366^^1 q^ + 0{q'^) . 



(B.l) 
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Following the procedure described in the main text, we can rewrite in terms of (quasi) 
modular functions according to 



3456 
1 

~ 12 



(4i? -s'^+p) (4i? - + 3p) 80R^ - 8R{5s^ - 22p) + bs'^ - Us^p + 9 V 
80R^ - 8i?(5s^ - 22p) + bs"^ - AAs'^p + 9^1 (Ti^^ - Ta^f) 



El 



+ {4R - + 6p) 64R^ - 80^2 (35^ - 4p) + 4i?(27s'^ - SSs^p + 49/ 



- 13s^ + 68s^p - 94s V + 30p^ + 2304iV 



E2 E4 



-^{4R-s^ + 6p) {R-s^+ p) \Tiel{2ei + Ot) + Ts^f (0^ + 2^1) 



E2 



+ 8{Tiel - T2e^fE2 + 16 Tt{el + 2ef,el - Ti{et + 2el)el - T1T2 etei{et - el) 
1 



+ 



2816i?^ - 30464i?^s^ + 3M24R^p + 67872^^5^ - 190848i?^s^p 



16i?^ - 8R{bs' - 6p) + 21s^ - ms^p + 3 V {Tiet - Ts^s) Ea 
1 

120960 

+ 94304i?V - 28400i?s^ + 134112i?s''p - 165664ii:s2p2 + 47616i?/ 
- 774144Ar(s2 -p) + 331776NR - 552960(rf + T1T2 + T^) 
+ 3323s^ - 22216s^p + 46862sV - 34584s V + 6615p''l Eq . 



(B.2) 



By expanding the modular functions in powers of g as shown in Appendix A, one can 
recover the instanton terms in (B.l). 



C. Reformulating the moduleir anomaly equations 

The modular anomaly equation (3.30) implies the following relation: 

n 1 2 , ^162 



(C.l) 



where the functions <^^'s have been defined in (4.13). Since the operators D and {—da) 
commute, by applying the latter to (C.l) and remembering that (pt+i = —da(p£, it is 
straightforward to obtain the action of D on any cpg and verify that these functions form a 
ring under it. For example, at the next step we obtain 

D(pi = (pi(p2 + ^ (P3 = (pi (- da(pi) + ^ dl(pi . (C.2) 

which is, up to rescalings, the viscous Burgers equation 

dtu = udxU + vd^u (C3) 

with the viscosity u proportional to eie2- 

In Section 4 the S-duality requirements were formulated in terms of the quantities S^^^ 
defined in (4.23), which we rewrite here for convenience as follows: 



00 ^ 



n=0 



(C.4) 
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with 

C=<pi 

Prom this definition it is easy to show that the following relation holds: 
In turn, this implies that 

27riro 



(C.5) 



(C.6) 



(C.7) 



Recalling that 27riro — (^2 = 27riT (see (4.14c)), this is tantamount to the recursion relation 
(4.24), namely 

(C.8) 



which allows to easily obtain the expression of any S^^^ once E^^^ is known. 

To determine E*^") we start considering the quantities Pi^\ Clearly, from (C.5) we 



have Pq^^ = ipo. The next case is 

= (D + caa)(^o 



n 2 1 2 , ^1^2 

C=V5l ^ ^ 



(C.9) 



where the last step follows from (C.l). By further applications of the operator (-D + C dg) , 
we get 

^2 = ^ ¥^2 + ( ^ ) V>4 , 

2^2, , (C.IO) 



eie2</?2 + (— j (6V329'4 + 5(^3) + j (/?6 , 

with similar expressions for higher values of n. In fact, it is possible to derive the general 
expression of the Pn^^^s at the first orders in their explicit dependence on eie2'- 



p(0) 
n 

p(0) 
n 



(n-1)! 



2 ' 



(n-1)! /I / 9 \2 



(C.ll) 



2^H9^2 



Inserting these results into (C.4), we obtain 

S(0) 

e(0) 



^° 2(27riro) 



= — log 1 
1 2 



1 ^4 



+ 



27riro 
5 



- log — 
2 ^ro 



(C.12) 



2 8(27rir)2 24 (27rir)3 ' 
where from the second equality on we have used that 2mTo — (p2 = ^ttit. 
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Applying the recursion formula (C.8) we can easily get the explicit expressions for the 
E^^^'s that are needed in the calculations presented in Section 4 or in Appendix D. They 
are, for £ = 1: 

1 'P3 



2 (27rir)2. 

l (P5 2 ifsipi 5 



^'''^°(8(27rir)3 + 3(27riT)4 ' 8 (27riT) 



+ 



for e = 2: 



S(2) 



27riTo 



E(^)|^ = (2.iro)^(i 



^2 

27rir / 

1 (^4 



(27riT)3 (27rir)4/ ' 



for £ = 3: 



for ^ = 4 



^^'''^°^'(2 (27rir)4 ' 2 (27rir)5 ' ^ (27rir) 
S(^)[ = (27riro)^( 



+3. '^3 



and finally for £ = 5 : 



S(5)|^ = (27rirof(- 



.(27riT)4 (27rir)5 



+10,^ + 15 



(27riT 



. (27rir)5 (27rir)6 

D. S-duality at order (£162)^ 

We start from the relation (4.18b) for n = 2, namely 



ToS{X2) ^-X2 = -ToS{ipi) ^ - ToS{Xi) 



^1 

irlV 



From (4.22) we read that 



To 5(991 



2 27riTn 



X2 

1 2mTn 



+ ^ 



1 ^1 5^(3) 



2 (27riro)2" 



(CIS) 



(C.14) 



(C.15) 



(C.16) 



(C.17) 



(D.l) 



(D.2) 



Substituting the expressions of the various 's given in Appendix C, and that of Xi 

k 

given in (4.31), we get 



ToS{ipi) ^ = 27riro(^- 



1 (p5 13 ip4:(p3 13 (fil \ _ _if2_ 

(27rir)3 24 (27rir)4 32 {2mT)^ ) 2mT 



X- 



2 ■ 



On the other hand, from (4.31) it follows that 



^'^^^'^-Sniri Sir) 



(D.3) 



(D.4) 
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The numerator of this expression can be computed from (4.22), yielding 



ZTTlTo 



+ eie2 



1 27riro 



+ 



(27riro)=^ j^;^ + €162 



(27riro)^(- 



1 (P5 13 V93V?4 , 13 (fl 



(27rir)3 ' --^-^ v—^'u/ ^2 (27rir)4 ^ 4 (27rir)5 4 (27rir) 



(D.5) 



+ 



where the last step follows from the results given in Appendix C. For the denominator we 
take into account that 



5(r)=5(ro)- 



1_ S{ip2) 
To 



27ri To 27ri 

Resorting again to (4.22) to evaluate Tq 5(922); we obtain 



(D.6) 



5(r) = -- 



1 



To 27riro 



S(2) 



.S(3) ) 
0/ 



/I 



---6162^- 



3 

+ T 



1 27riro o> 

)..... 



+ 



(D.7) 



.2 (27ri)2T3 4 (27ri)3r4 
Inserting (D.5) and (D.7) into (D.4) and extracting the term of order 6162, we find 

ro5(Xi) ^ = -(27nro)(- + - + - . (D.8) 

Using this result and (D.3) into (D.l), we finally obtain the following constraint on X2: 



7 (/33V34 7 



(D.9) 



Considering the structures involved in this relation, we can try an Ansatz such that X2 is 
written as a total derivative: 



X2 = da( -Ai 



(27rir) 



-A, 



(27rir)= 



A,,^ + 2(Ai + A2),^+3A2: 



(D.IO) 



(27rir)2 ^ ^ ^'(27rir)3 ^ (27rir)4 ' 
With this position we find straightforwardly, using the formulae of Appendix C, that 



r,S(X,)l - 2.ir„(^) = 2.ir„((6A, - 4A.) ^ + (9A, - 6A,) ^) . (D.ll) 

Comparing this result to the constraint (D.9), we obtain a single independent equation: 

6Ai-4A2 = ^. (D.12) 

Our Ansatz thus satisfies the requirement (D.9) with still one free parameter. In Section 5 
we chose to fix this arbitrariness by requiring that the modified prepotential F introduced 
in (5.2) differs from F by a series of corrections in powers of €162 whose coefficients only 
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involve (/33 and the modified effective coupling r defined in (5.5). It is easy to see that, 
starting with X2 given in (D.IO), the F term of order (6162)^ in (5.7) would contain a part 
proportional to (^4 unless Ai assumes a specific value which, by virtue of (D.12), fixes also 
A2. These values are: 

which are those used in (5.1). 
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